Decoherence and relaxation of a qubit coupled to an Ohmic bath directly and via an 
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Abstract 

Using the numerical path integral method we investigate the decoherence and relaxation of qubits coupled to an Ohmic 
bath directly and via an intermediate harmonic oscillator (IHO). Here, we suppose the oscillation frequencies of the bath 
modes are higher than the IHO's. When we choose suitable parameters the qubits in the two models may have almost same 
decoherence and relaxation times. However, the decoherence and relaxation times of the qubit in the qubit-IHO-bath model 
can be modulated through changing the coupling coefficients of the qubit-IHO and IHO-bath and the oscillation frequency of 
the IHO. 



PACS numbers: 67.57.Lm, 03.65.Yz, 31.15.Kb 

I. INTRODUCTION 

Since Shor's algorithm [l[ for factoring large numbers, 
the theory of quantum computation and quantum infor- 
mation has attracted great recent interest. It is believed 
that quantum computers may perform some useful tasks 
more efficiently than their classical counterparts. Despite 
the great promises of performing quantum computations, 
however, there are still many practical difficulties to be 
resolved before quantum computers might become avail- 
able in future. One of the difficulties is that the qubit 
has too short decoherence time, which is in fact a central 
impediment for practical qubit to be taken as the cell 
of quantum computers. Solid-state qubits are considered 
to be promising candidates for realizing building blocks 
of quantum computers because of their integrability and 
flexibility in the design. However, the quantum logical 
gates or registers made up with these qubits have still 
shorter decoherence and relaxation times. These mo- 
tivated a lot of studies on the decoherence, relaxation 
and manipulation of the qubits. Many significant results 
in the field, not only theoretical but also experimental, 
have been achieved. Most of the theoretical research is 
based on a spin-boson (SB) model which is supposed to 
be constructed with a spin (qubit or two-level system) 
coupled to a bath. The model has many physical cor- 
respondences and has been widely investigated in recent 
years 0, H, 0| ■ In order to investigate a qubit in its envi- 
ronment, there is another model which is different from 
the original spin-boson one. In the model, the qubit cou- 
ples to the coordinate X of a harmonic oscillator, which 
we shall sometimes call the "intermediate harmonic os- 
cillator" (IHO), and which in turn is coupled to a bath. 
We call the model the spin-IHO-bath (SIB) model. Re- 
cently, this model receives much interest in the context 
of quantum computing with condensed matter systems, 
especially with superconducting flux qubit devices, see 
Ref. [|[ and within. It is also useful to investigate the 
measurement of the solid-state qubit Q and magnetic 
resonance force microscopy @, [U • 

On the other hand, the SB and SIB models can be used 



to describe the electron transfer in chemical and biolog- 
ical molecules. Many theoretical investigations in this 
field use the SB model @, 0, El El ■ Based on Leggett 
[lH, Garg et al. investigated the SIB model, and ob- 
tained a map of it to the SB one, and obtained the spec- 
tral density of the effective bath for the map. They used 
this model to study the migration of an electron from one 
biomolecule to another, or between two localized sites in 
the same biomolecule. It has been shown that the co- 
herence is very important not only to qubits for making 
quantum computers but also to electrons for transferring 
energy in biological systems flH [T3 |. However, if one 
wants to know the decoherence and relaxation behaviors 
of the two-level systems in their environment, essentially, 
the dynamics of the systems needs to be known. 

If the qubit energy splitting (denoted by A hereinafter) 
is not equal to zero, the two models are not exactly solv- 
able. However, they can be analyzed using adiabatic 
renormalization in which a systematic weak damping ap- 
proximation must be used. They can also be investigated 
with some approximation methods based on the pertur- 
bative scheme which also asks for the systems (qubits) 
weakly cou plin g to their environment. Many other meth- 
ods llTl. [TsL Tl9l] for solving the models have been proposed 
in recent years, most of which are based on the Born- 
Markov approximation. However, it has been pointed 
out that the use of the approximation is inappropriate 
at the large tunneling amplitude and low temperatures. 
Recently, some different schemes to solve the SIB model 
have also been put forward. Gassmann et al. 8j obtain 
an exactly solvable model from the SIB through dropping 
an unimportant term, where an approximation similar to 
the Born-Oppenheimer one is used. This method is not 
successful as the "dropping term" is not small enough. So 
it is important to find out some methods to accurately 
estimate the dynamics of the qubits in the two models. 
Based on the insight into the dynamics we may under- 
stand the decoherence and relaxation better and may 
bring forward some schemes on how to suppress them. It 
is also of much interested to find out the qubit in which 
model, SB or SIB, has longer decoherence and relaxation 
times. These problems are not obvious. 
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An excellent method, accurate numerical path integral 
method based on the qusiadiabatic propagator path in- 
tegral (QUAPI) scheme [13, [H[ may be a suitable tool 
for solving the two models. To our problems we choose 
the iterative tensor multiplication (ITM) algorithm for 
the numerical scheme. As Makri pQ , l2ll ] addressed that 
the method is non-Markovian and it can make the calcu- 
lations accurate enough even at very low temperatures, 
large tunneling amplitude and strong couplings for which 
the Markovian approximation is unsuitable. In this pa- 
per we shall use the tool to investigate the dynamics and 
then the decohcrcnce and relaxation of the qubit in the 
SB and SIB models. 



II. MODELS AND DYNAMICS 

The Hamiltonian of the SB model is 

2 1 



e A 
Hsb = h I + 2" Gx 



-E 



Pi 

2m, 



x t + 



Suppose the bath has an Ohmic spectral density 

r2 tt 



J ohm, M = 77 E — & ( W 

2 i rriiUJi 



— Uj/uj c 



mnoj 
(1) 



(2) 



Here, £ is the dimensionless Kondo parameter [12|, [22| 
(the relationship of £ with the friction coefficient r\ is 
£ = 2ri/nh 23]), <Tj (i = x, z) are the Pauli matrix, and 
to c is the high-frequency cutoff of the bath modes. This is 
a well-known quantum dissipation model and it has been 
widely investigated 0, Q • 

If we consider the qubit coupling to the coordinate X 
of a single IHO which in turn is coupled to a bath, and 
if we let the couplings be linear, the Hamiltonian of the 
SIB system reads 
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where M and P are the mass and momentum of the IHO, 
and the displacement A characterizes the coupling of the 
qubit to the IHO, and KCj are the coupling coefficients of 
the IHO to the bath modes. It is shown that the system 
has a one to one map to the following system [lj| 
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with an effective spectral density (see Appendix) 
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where V = k 2 t]/2M. When the bath modes have lower 
frequencies than CIq, the dynamics of the qubit in the 
bath with effective spectral density J e ff (w) is similar to 
the dynamics of the qubit in Ohmic bath with spectral 
densit y J n hm (w) . which is widely investigated in recent 
years [H HE Hfl [13, H S Hf. In this paper we in- 
vestigate another limited case, namely, the bath modes 
have higher frequencies than the oscillation frequency Hq 
of IHO. The length of the memory times of the baths can 
be estimated by the following bath response function 
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Here, (3 = l/fc^T where ks is the Boltzmann constant, 
and T is the temperature. It is shown that when the real 
and imaginary parts behave as the delta function 5 (t) 
and its derivative 6' (t) , the dynamics of the reduced 
2 density matrix is Markovian. However, if the real and 
imaginary parts are broader than the delta function, the 
dynamics is non-Markovian. The broader the Re[a (t)] 
and Im[a (i)] are, the longer the memory time will be. 
The broader the Re[a (t)] and Im[a (t)} are, the more 
serious the practical dynamics will be distorted by the 
Markovian approximation. The memory time of the ef- 
fective bath is determined by T. The larger the T is, the 
shorter the memory time of the effective bath will be. 
If the memory time of the bath is short enough, the re- 
duced density matrix may be obtained in virtue of the 
Markovian approximation. However, if the memory time 
of the bath is too long, the ITM method is in fact in- 
appropriate, saying nothing of other methods based on 
Markovian approximation. Clearly, the value of the T 
may vary according to the difference of the physical sys- 
tems. For example, when the persistent-current qubit is 
measured by a dc SQUID, the system can be modeled 
by Eq. © with Eq. ©, here T = 1/R S C S . Typically, 
R s = 100 n, C s = 5 pF, so T - 10 u , see Ref. In this 
paper we set T = 2.6 x 10 11 . This value of T makes the 
calculation with ITM suitable but other methods based 
on Markovian approximation unsuitable. 

In Fig. 1 we plot the Re[a (t)] and Im[a (t)] as J (to) = 
Johm(^), and J e ff (to), where we set Xk = 1050, £ 
0.01 . n = 10A, T = 0.01, r = 2.6 x 10 u , A = 5 x 10 9 
Hz, the lower-frequency and high-frequency cut-off of the 
baths modes u>q = 11A, and to c = 100A. It is shown that 
the memory time of the Ohmic bath is similar to the one 
of the effective bath, namely, lo c t™b ~ ^c^sib- Both 
baths have shorter effective memory times. This will be 
rificd in Fig. 2 below. 

The dynamics of the qubit is characterized by the time 
evolution of the reduced density matrix, obtained after 
tracing out the bath degrees of freedom, i.e., 

p (s", s'; t) = Tv bath (s"\ e- mt ' h R (0) J m ' h \s') . (7) 

In actual cases the initial state of the total system must 
be entangled even at the beginning of the evolution. 
However, at the beginning the entanglement is very weak, 
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otherwise the qubit must have lost its coherence and re- 
laxed to its thermal equilibrium in a very short time. 
This is not our case investigated in this paper. So, for 
simplicity of the presentation we assume that the entan- 
glement (through interaction) between the qubit and its 
environment is switched on at t = 0, i.e.. the initial den- 
sity matrix has the form [13, [13, S lH 



R(0) = p(0)® Pbath (0), 



(8) 



Y a x) and H env = Hsb — Hq, 
In order to make the calcula- 



where p (0) and p 6ath (0) are the initial states of the 
qubit and bath. The scheme which we use to calcu- 
late the reduced density matrix p(t) is a well established 
ITM algorithm derived from the QUAPI. It is a nu- 
merically exact algorithm and is successfully tested and 
ado pted in various problems of open quantum systems 
fl2l |22L I23I ]. For details of the scheme, we refer to pre- 
vious works [13, Hi]]. The QUAPI asks for the system 
Hamiltonian splitting into two parts Hq and H env . Here, 
we take H = H (^<J Z H 
or H env = Hsib — Hq 
tions converge we use the time step uj c At = 0.6, which is 
smaller than the characteristic times of the qubits in the 
systems. 



III. DECOHERENCE AND RELAXATION 

The decoherence is in general produced due to the 
interaction of the quantum system with other systems 
with a large number of degrees of freedom, for exam- 
ple the devices of the measurement or environment. To 
measure the decoherence one may use the entropy, the 
first entropy, and many other measures, such as the 
maximal deviation norm, etc. (see for example Refs. 
[3ll I32I HH ). However, essentially, the decoherence of 
an open quantum system is reflected through the decays 
of the off-diagonal coherent terms of its reduced den- 
sity matrix [34j |. The decoherence time denoted by t 2 
measures the time of the initial coherent terms to their 
1/e times, namely, p i (n,m) p* (n, to) = Pi (n, to) /e. 
Here, n ^ m, and n, to = or 1 for qubits. In the 
following, we investigate the decoherence via directly 
describing the evolutions of the off-diagonal coherent 
terms, instead of using any measure of the decoherence. 
Similar to the decoherence, the relaxation of the qubit 
can also be investigated with the diagonal elements of 
the reduced density matrix. The relaxation time is de- 
noted by T\, which measures the time of an initial state 
to the final thermal equilibrium state through estimat- 
ing the diagonal terms of the reduced density matrix, 



namely, p i (n, n) 
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In the following calcu- 



lations we assume the initial state of the environment 
Pbath (0) = nfc e_/3Mfc / Tr fe (e-^ M ") . As calculating the 
off-diagonal element p 12 we let e = 10A which can make 
the p 12 decay stably. If e — ■+ A the p 12 will decay with 
some oscillations, which may affect our judgement on de- 
coherence times. The closer the two parameters are, the 
more strongly the matrix elements will oscillate. When 



we calculate p lx we choose parameters e = A because 
the oscillations of the p u do not affect our judgement 
on relaxation times from the figures. The reader should 
note that, the increase of the A and e will shorten the de- 
coherence and relaxation times, so the decoherence time 
t 2 and the relaxation time t\ in the figures are not com- 
parable because they are plotted in different two sets of 
the parameters. But, it is clear that the relaxation time 
T\ is longer than the decoherence time t%. 

In the ITM scheme, one should at first choose the fc max , 
so that fc max w c Ai must be larger than the effective mem- 
ory time u} c t™ b and w c r™ JB of the baths. In Fig. 2 we 
plot the reduced density matrix elements p n and p 12 as 
femax = 2, 3, and 4 for the SB and SIB models. Here, we 
set Xk — 1050, the low- and high-frequency cut-off are 
ujq = HA, and oj c = 100A. From the Fig. 2 we obtain 
that as Xk — 1050 the qubit has almost the same deco- 
herence and relaxation times in SB and SIB models, and 
the calculations are in fact convergent as fc max ^3. It 
is known that the qubit will show different decoherence 
and relaxation when it has different initial states. We 
can not obtain general results on the decoherence and 
relaxation from some special initial states. Essentially, 
this limitation is inevitable, because in order to obtain 
measures of the decoherence and relaxation one should 
perform the calculations over all possible initial states in 
general [30l [3l| . But, we may figure out some informa- 
tion on the trend of the decoherence and relaxation of 
the qubit from some special results. In particular, we 
are more interested in some special results because some 
schemes on quantum information are based on the usage 
of some special initial states. In Fig. 3 we plot the de- 
coherence and the relaxation of the qubit in SB (a) and 
SIB (b) models in different initial states. These states 

/ l72, 

|C 7 ) = [1, 0] T . As Xk > 1050 (or Xk < 1050) the decoher- 
ence and relaxation times will be shortened (or length- 
ened), as shown in Fig. 4. The frequency of the IHO will 
strongly affect the qubit decoherence and relaxation. If 
the Slo of the IHO decreases, the decoherence and relax- 
ation times will be lengthened, and vice versa. Namely, 
the farther the f2o departs from the lower-frequency cut- 
off of the bath, the longer decoherence and relaxation 
times the qubit will have, as evidenced in Fig. 5. 



IV. CONCLUSIONS 

In this Letter we have investigated the decoherence and 
relaxation of a qubit coupled to an Ohmic bath directly 
and via an IHO. In our investigations, we fix the tun- 
neling splitting A of the qubit and assume that the IHO 
is far off the resonance to the environment modes, i.e. 
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the oscillation frequency of the IHO is smaller than the 
lower-frequency cut-off of the bath modes. In addition, 
we suppose that the damping coefficient T of a bath to the 
IHO is moderate. If the quantity is too small, the mem- 
ory time of the bath may be very long, thus the methods 
based on not only the Markovian approximation but also 
the ITM scheme are inappropriate for the investigation of 
the qubit dynamics. On the contrary, if the T is too large, 
and the memory effects of the bath can be neglected, the 
dynamics of the qubit can be investigated in many other 
ways based on the Markovian approximation. In spite 
of these limitation our investigations have wider appeal 
because the method is suitable for a wider range of T, 
and most practical physical systems have the moderate 
r. By use of the accurate numerical path integral scheme, 
the ITM based on the QUAPI method we have obtained 
the evolutions of the reduced density matrix elements of 
qubit in the SB and SIB models. We have obtained that: 
(1) When Xk = 1050 the qubit in the SIB model has 
almost the same decoherence and relaxation times with 
the qubit in the SB model. If Xk > 1050 the damping 
of a bath to the IHO or (and) the IHO to the qubit in- 
crease, which results in the decrease of the decoherence 
and relaxation times, and vice versa. (2) The decoher- 
ence and relaxation times of the qubit in the SIB model 
increase with the decrease of the oscillation frequency ilo 
of the IHO, and vice versa. (3) The decoherence and re- 
laxation times of the qubit in the SB and SIB models will 
increase with the decrease of the e and A, which has not 
been plotted in the Letter. The longer decoherence and 
relaxation times are necessary for not only the qubits for 
making the quantum computers but also the electrons 
for transferring energy in biological systems. In order 
to make the qubits or electrons in the SIB model have 
longer decoherence and relaxation times we may try to 
make the Hq an d Xk smaller. 
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V. APPENDIX 

In the appendix we derive the effective spectral density 
expressed with Eq. ©. As Ref. [lU pointed out that 
the map from Eq. [3] to Eq. [4] does not involve the spin 
(qubit), the same J e // ( w ) w ih control the dynamics of 
a continuous variable q moving in some potential U(q) 
and coupled to coordinates X and {xt} in the same way 
as the spin. So we can deduce the J e ff(ui) with the 
Hamiltonian 



H = 
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where p q is the momentum conjugate to q. Defining U' = 
dU jdt and using the dots for the time derivatives, the 
classical equations of the motion are 
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Using the Fourier transforms, Eqs. (I10M12I) can be writ- 
ten as 
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Insetting Eq. (fT~5|) into Eq. lfT4|) . we have 
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and Eq. (|16j) we have 
X(lu) 
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Insetting Eq. (Ql]) into Eq. ((131), we nave 

r MVLqX 2 [Mil 2 + L(lu)] (Mil 2 ,) 2 X 2 
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= -K(q). 
Introducing the function K(u), and making 

mh 2 x 2 l(lu) 
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By using Eqs. ([2]) and (fT6|) we have 
L (uj) = -lu 2 ( M + k 2 



q(Lu) 
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Taking the cut-off frequency to be infinity, we have Here, T = K 2 r)/2M, and £ = 2rj/irh. 

L (u) = -Muj 2 + ir]K 2 co. (22) 

Substituting this in Eq. ([20]) . from J e ff (u) — 
lim Im [K(u> — is)] (u is real), we have 

Jeff H = |aV^ - "| (23) 

2 (w 2 - fig) + 4T 2 lo 2 
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VI. FIGURE CAPTIONS 

Fig. 1 The response functions of the Ohmic bath and 
effective bath, where A = 5 x 10 9 Hz, Xn = 1050, £ = 
0.01, tt = 10A, T = 0.01, r = 2.6 x 10 11 , the lower- 
frequency and high-frequency cut-off of the baths modes 
luq = 11A, and ui c = 100A. 

Fig. 2 The evolutions of reduced density matrix el- 
ements p 12 and p lx in SB and SIB models in different 
values of K max . Here, e = 10A for p 12 , e — A for p n , 
the initial state is 1^} and the other parameters are the 
same as in Fig. 1. 

Fig. 3 The evolutions of reduced density matrix ele- 
ments of p 12 an d Pn in SB (a) and SIB (b) models in 
different initial states. These initial states are given in 
the context. Here, -Kmax = 3, and the other parameters 
are the same as in Fig. 2. 

Fig. 4 The evolutions of reduced density matrix ele- 
ments of p 12 and Pn in SIB model in different values of 
k\, the other parameters are the same as in Fig. 2. 

Fig. 5 The evolutions of reduced density matrix ele- 
ments of p 12 and Pn in SIB model in different values of 
f^O; the other parameters are the same as in Fig. 2. 
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